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A  zero-dimensional  time-dependent  high-pressure  (5-106  -  5-108  Pa)  slab  capillary 
discharge  model  is  presented.  The  model  includes  a  heat  transfer  radiation  model  based  on  a 
radiation  database.  This  database  has  been  constructed  using  commercially  available 
radiation  software  PrismSPEC  to  calculate  the  radiation  heat  flux  output  from  an  uniform 
plasma  slab.  Thus,  unlike  earlier  models,  this  model  does  not  use  any  “asymptotic”  radiation 
models,  but  self-consistently  calculates  the  radiation  heat  flux  at  the  thin  transition  layer, 
between  the  uniform  plasma  core  and  the  ablative  capillary  walls.  The  model  includes  the 
thermodynamics  of  partially  ionized  plasmas  and  non-ideal  effects  taking  place  in  the  high- 
density  plasma  and  assumes  local  thermodynamic  equilibrium  (LTE),  fully  dissociated 
plasma,  no  heat  losses  into  the  capillary  walls,  a  ratio  of  thermal  pressure  to  magnetic 
pressure  much  larger  than  unity  (P  »1),  and  the  existence  of  a  sonic  condition  at  the  exit 
plane  (the  plasma  flow  is  expected  to  be  chocked  at  the  bore  exit).  The  model  predicts  the 
existence  of  two  steady-state  regimes  of  plasma  pressure  for  ablative  discharge  operation  at 
a  given  plasma  temperature.  The  first  regime  occurs  when  the  plasma  is  so  dense  (~  1026  m'3) 
that  the  radiation  mean  free  path,  Amd  is  less  than  the  slab  gap  of  capillary,  Da,  the  case  of 
super-high  pressure  capillary  discharge.  The  second  regime  occurs  when  the  plasma  density 
is  much  lower  (~  1024  -  102S  m'3)  such  that  Arad  is  much  larger  than  the  capillary  gap,  i.e.  the 
case  of  moderately  high  plasma  pressure.  Both  regimes  converge  at  small  plasma 
temperature,  and  there  is  no  steady-state  solution  for  small  plasma  temperatures.  The 
calculations  show  that  with  an  increase  in  the  capillary  length,  the  density  of  the  plasma  in 
the  second  regime  (Lrad  »  Da )  may  become  so  small,  that  the  conduction  heat  flux  to  the  wall 
becomes  larger  than  the  radiation  flux  to  the  capillary  wall  and  transfers  from  the  radiation 
mode  to  thermal  conduction  mode.  We  show  that  in  this  mode  the  capillary  discharge  can 
still  exist  when  LTE  is  valid  and  p  »1.  In  the  Lrad  »  Da  regime  the  radiation  heating  of  the 
wall  for  all  capillary  lengths  is  always  a  few  orders  of  magnitude  larger  than  conduction 
heating  of  the  wall.  Both  regimes,  radiation  and  thermal  conduction,  may  be  attractive  for 
thruster  applications  depending  on  specific  configurations. 


I.  Introduction 

The  role  of  space-based  systems  for  both  commercial  and  government  customers  continues  to  evolve, 
continually  providing  new  requirements  for  the  development  of  satellite  propulsion  systems.  One  development  path 
is  towards  ever  increasing  power  levels  in  a  compact,  efficient  thruster.  Being  very  efficient,  high-pressure  ablative 
capillary  discharge  is  a  good  candidate  for  high-power  plasma  thrusters.  Creating  a  comprehensive  model  of 
capillary  discharge  is  important  to  understand  the  physics  and  engineering  aspects  of  the  capillary  discharge 
thruster. 

A  schematic  of  a  slab  capillary  discharge  thruster  is  shown  in  Fig.  1,  where  Dc  is  the  slab  gap  and  L,  is  the 
capillary  length.  The  discharge  maintains  a  resistive  arc  through  a  narrow  insulating  capillary  by  the  continual 
ablation  of  the  capillary  wall  material  as  shown  in  Fig.  1  or  by  injected  mass.  The  ablative  capillary  discharge  can  be 
described  via  a  three-layer  configuration.  The  outermost  layer  is  a  solid  wall,  usually  some  form  of  polyethylene, 
occasionally  Teflon  or  some  other  insulating  material.  Material  from  the  wall  is  evaporated  and  enters  the  thin 
transition,  or  “boundary”  layer,  where  it  is  dissociated,  ionized,  and  heated  to  plasma  temperature.  The  innermost 
layer  is  the  plasma  core.  The  closed  end  of  the  capillary  (left  side  of  Fig.  1)  is  one  electrode;  the  other  electrode 
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(right  side  of  Fig.  1)  is  an  open  end,  through  which  the  plasma  can  flow  and  expand.  The  polarity  has  no  appreciable 
effect  on  discharge  operation.  Ohmic  heating  is  responsible  for  heating,  ionization,  ablation,  and  radiation.  In  the 
model,  a  sonic  condition  ( Ma  =  1,  in  Fig.  1)  is  assumed  to  exist  at  the  open  end  of  capillary. 


In  terms  of  plasma  propulsion  concepts,  the  capillary  plasma  thruster  has  to  satisfy  two  main  conditions.  First, 
the  capillary  plasma  should  be  in  a  local  thermodynamic  equilibrium  (LTE),  indicating  that  all  sorts  of  particles, 
electrons,  ions  and  neutrals,  have  the  same  temperature.  Then,  the  electrical  current  will  heat  both  the  elections  and 
the  heavy  particles  (i.e.,  not  only  the  electrons),  thus  providing  high  thrust  per  watt.  Second,  as  it  is  well  known, 
pinching  the  capillary  discharge  leads  to  a  formation  of  a  narrow  plasma  core  region  with  high  electron  temperature 
and  a  large,  relatively  cold  peripheral  plasma  region.  This  also  decreases  the  efficiency  of  the  capillary  thruster 
because  only  a  small  amount  of  input  discharge  energy  is  transferred  to  the  heavy  particles,  with  almost  all  of  the 
input  energy  being  transferred  to  the  electrons.  To  minimize  the  pinching  of  the  capillary  discharge,  the  ratio  of  the 
thermal  pressure  to  the  magnetic  pressure,  /(.  in  the  capillary  plasma  thruster  should  be  greater  than  unity. 

The  ablative  capillary  discharges  that  have  been  previously  studied  as  potential  electro-thermal  and  electro- 
thermo-chemical  gun  devices  by  various  investigators  since  the  mid-1980s  see  111  and  references  therein,  satisfy 
both  of  these  conditions  well.  These  studies  investigated  the  dynamics  of  high-pressure  ablative  plasma  discharges; 
they  gave  pressures  in  a  range  of  0.1-1  GPa  and  temperatures  of  the  order  of  1-3  eV.  Only  in  recent  theoretical 
papers1”"14  have  the  authors  explored  applications  of  such  discharges  to  electro-thermal  capillary  thrusters. 

In  all  theoretical  models  of  ablative  capillary  discharges1"7' 1015  zero-  and  one-dimensional,  the  thickness  of  the 
transition  layer,  Fig.  1,  is  assumed  to  be  negligibly  small  compared  to  the  inner  capillary  radius  and  the  radial 
distributions  of  plasma  parameters  in  the  plasma  core,  Fig.  1,  are  assumed  to  be  uniform.  This  statement  is  based  on 
the  premise  that  the  radial  radiation  heat  transfer  across  the  capillary  chamber  is  so  strong  that  it  flattens  the 
temperature  and  density  distributions  across  the  plasma  core  leading  to  a  very  thin  transition  layer.  As  it  has  been 
mentioned  in  2,  a  strong  plasma  convection  that  might  develop  in  the  capillary  discharge  can  also  lead  to  flattening 
the  capillary  discharge  parameters  in  the  plasma  core  and  thin  boundary  layer.  Unfortunately,  there  exists  no  model 
of  the  ablative  capillary  discharge  that  would  verify  this  premise.  The  main  challenge  for  constructing  such  a  model 
is  a  comprehensive  description  of  the  radiative  heat  transfer  in  the  capillary  discharge. 

It  is  worth  noting  that  in  other  types  of  capillary  discharges  operating  in  non-ablation  regime  with  much  faster 
pulses,  lower  pressure,  high  electron  temperature  in  the  plasma  cores  (no  LTE)  and  /?  <  1 .  the  plasma  density  is 
usually  so  small  that  thermal  conduction,  not  radiation,  is  responsible  for  the  heat  transfer  across  the  capillary  and 
the  creation  non-flat  radial  plasma  density  and  temperature  distributions16"19.  These  discharges  have  been  used  as  X- 
ray  radiation  sources  in  spectroscopy  and  lithography  for  wakefield  acceleration  and  ultra-intense  laser  guiding. 
Unlike  in  ablative  capillary  discharges,  in  these  plasma  sources,  capillary  chambers  are  filled  in  by  a  gas  and  the 
capillary  walls  are  made  from  a  ceramic  material  preventing  large  ablation.  However,  as  it  will  be  shown  in  this 
work,  ablative  capillary  discharges  may  also  be  operated  in  low  pressure  regime  with  thermal  conduction  heat  larger 
than  radiative  heat  (as  in  non-ablative  capillary  plasma  sources),  but  with  plasma  in  LTE  and  /?  >  1  as  desired  for 
propulsion  applications. 

The  first  comprehensive  zero-dimensional  analytical  model  of  ablative  capillary  discharges  was  developed  by 
Loeb  and  Kaplan2.  By  balancing  Ohmic  heating  by  the  radiation  losses,  wall  ablation  and  plasma  outflow,  they 
obtained  scaling  laws  for  the  plasma  temperature,  pressure,  wall  ablation  rate,  and  other  parameters  of  the  capillary 
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discharge  as  functions  of  the  current,  the  capillary  radius  and  length.  Some  key  assumptions  made  in  this  model 
include:  (1)  plasma  radiation  is  modeled  as  blackbody  radiation;  (2)  the  radiation  flux  incoming  to  the  boundary 
layer  (Fig.  1)  ablates  the  wall  material,  dissociates,  ionizes  and  heats  ablated  vapor  to  the  plasma  conditions,  and  no 
heat  losses  occur  in  the  bulk  wall;  (3)  the  plasma  is  fully  ionized;  (4)  when  calculating  the  plasma  enthalpy,  the 
authors  neglect  the  energy  cost  of  ionization.  The  last  assumption  can  be  a  serious  error,  since,  for  example,  the 
carbon  ions  can  be  twice  and  even  three  times  ionized  and,  therefore,  the  contribution  of  the  ionization  cost  into  the 
total  plasma  enthalpy  can  be  very  large,  causing  the  dependence  of  enthalpy  on  temperature  to  be  very  different  than 
assumed  in  model2.  These  simplifications  allowed  Loeb  and  Kaplan  to  obtain  analytical  solutions  to  the  scaling 
laws,  a  significant  progress  at  the  time.  However,  these  scaling  laws  are  somewhat  questionable,  due  to  the  level  of 
approximations  made.  In  further  work3"7'  10"16,  the  Saha  equation  is  used  to  more  accurately  calculate  the  plasma 
composition  and  enthalpy  of  the  plasma. 

In  the  time-dependent  Calligan-Mohanti  zero-dimensional  model',  the  rate  of  mass  ablation  M  is  determined 
by  converting  the  radiation  incident  on  the  solid  wall: 

M-Kap=FraJw-Aw  ,  (1.1) 

where  hvap  is  the  specific  heat  evaporation  of  the  wall  material;  Aw  is  the  surface  area  of  the  capillary  tube;  Fradw  is 
radiation  heat  flux  incident  on  the  insulator  wall.  Gilligan  and  Mohanti  assume  Fmdw  to  be  a  fraction  f  (the  grey 
factor)  of  the  blackbody  radiation  leaving  the  arc: 

Frad,w=f-°sB-T*  ,  (1.2) 

where  <rSB  is  the  Stefan-Boltzmann  constant  and  T  is  the  temperature  of  plasma  core,  Fig.  1.  Assuming  local 
thermodynamic  equilibrium  and  using  the  Saha  equation,  Gilligan  and  Mohanti  calculate  the  internal  energy  of  the 
plasma,  the  resistivity,  and  other  thermodynamic  parameters.  To  get  agreement  with  experiments  they  must  assume 
that  only  a  small  fraction  (f «  1)  of  the  plasma  radiation  reaches  the  wall.  Thus,  this  model  demonstrates  the 
obvious  fact  that  only  a  small  fraction  of  the  energy  that  has  been  input  into  the  discharge  is  actually  spent  for 
evaporation  of  the  wall  and  the  main  part  of  input  energy  is  spent  on  dissociating,  ionizing  and  heating  the  ablated 
material  to  plasma  temperature. 

In  one-dimensional  models4"8  1011  and  references  therein,  the  authors  calculate  the  plasma  composition, 
resistivity,  temperature,  and  other  plasma  parameters  along  the  capillary  discharge  by  using  different  levels  of 
approximations  and  simplifications.  However,  all  these  models  assume  LTE  in  the  plasma  layer  and  »  1  and  no 
heat  losses  in  the  bulk  of  the  capillary  walls.  In  4"7'  10-1 1  the  authors  calculate  the  ablation  mass  density  rate  p 
averaged  over  the  cross-section  of  the  capillary  as 


P  = 


2 -f-aSB-T4/Rc 


£  +  P/  p 


(1.3) 


where  s,  P.  p,  and  T  are  the  internal  energy,  plasma  pressure,  density  and  temperature  in  the  plasma  core  layer; 
fcTsB-T4  is  the  radiation  heat  flux  incoming  to  the  transition-boundary  layer;  and  R,  is  the  capillary  radius,  while 
Raja-Varghese -Wilson8  have  used  basically  Eqs.  (1.1)  and  (1.2). 

In  these  models  the  grey  factor  ranges  from  0.6-0.85  and  0.05-0.1  to  fit  experimental  data  depending  on  where 
authors  consider  the  radiation  flux  either  at  the  plasma-transition  interface4"7'  10-1  ,  Eq.  (1.3),  or  at  the  ablative 
surface8,  Eqs.  (1.1)  and  (1.2).  This  indicates  that  plasma  radiation  in  experiments4' 6  8  cannot  be  described  by  a  pure 
blackbody  radiation  approximation,  and  the  Rosseland  radiation  mean  free  path  ).rad  is  larger  than  the  capillary 
radius. 

In  works14,  15  the  authors  used  a  different  (not-blackbody-type)  radiation  model,  a  volumetric  bremsstrahlang 
(free-free)  radiation  model.  In  their  radiation  models,  they  completely  ignore  the  recombination  (free-bound)  and  in¬ 
line  (bound-bound)  radiation  processes.  However,  our  calculation  of  radiation  balance  using  PrismSPEC  radiation 
software20  shows  that  the  free-bound  and  bound-bound  electron  transitions  give  the  main  contributions  into  the  total 
radiation  balance  in  high-pressure  capillary.  Thus,  using  only  bremsstrahlang  radiation  in  modeling  high-pressure 
capillary  discharge  is  questionable. 
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In  a  zero-dimensional  ablative  capillary  discharge  model1"  the  authors  have  shown  that  steady-state  ablative 
capillary  discharge  in  fact  may  exist  in  two  regimes.  The  first  regime  is  called  a  super-high  pressure  (SHP)  regime, 
where  plasma  density  is  assumed  to  be  extremely  high  so  that  the  plasma  is  opaque  to  its  own  radiation;  in  this  case, 
a  blackbody  radiation  approximation  is  applied.  This  regime  corresponds  to  the  previously  developed  models  with 
the  grey  factor  equal  to  one.  In  the  second  regime,  the  moderately  high  pressure  (MHP)  regime,  the  plasma  density 
in  the  plasma  core  is  assumed  to  be  much  smaller  so  that  the  plasma  is  completely  transparent  to  its  own  radiation, 
but  the  transition  layer  is  assumed  to  be  opaque  to  the  plasma  radiation  and  to  absorb  all  incoming  radiation.  In  this 
regime,  the  radiation  model  has  included  bremsstrahlang  (free-free)  radiation  and  recombination  (free-bound) 
radiation;  the  bound-bound  radiation  has  been  ignored.  It  should  be  stressed  that  the  last  assumption  limits  the 
applicability  of  this  regime  to  the  case  of  partially  ionized  plasma  and  multi-ionized  fully  ionized  plasma,  where  the 
line  radiation  (bound-bound  radiation)  can  be  much  larger  than  continuum  radiation.  However,  the  results  obtained 
in  this  moderately  high  pressure  regime  are,  probably,  a  valid  guideline  for  the  case  of  pure  hydrogen  (a  propellant 
of  choice  for  high-performance  propulsion),  where  hydrogen  plasma  is  assumed  to  be  fully  ionized  and  bound- 
bound  radiation  is  absent.  The  prospects  of  ablative  capillary  discharges  for  propulsive  applications  are  also 
discussed  in  this  paper. 

Thus,  all  previously  developed  models  of  ablative  capillary  discharges  either  consider  the  asymptotic  radiation 
cases,  blackbody  and  volumetric  radiation  regimes  or  use  blackbody  radiation  with  a  grey  factor  to  obtain  agreement 
with  the  experiments;  the  experiments,  probably,  correspond  to  intermediate  regimes,  between  the  blackbody  and 
volumetric  regimes. 

In  a  recent  zero-dimensional  high-pressure  (5T06  -  5T0S  Pa)  capillary  discharge  model1'  the  authors  used  a 
radiation  heat  transfer  model  based  on  a  radiation  database.  This  database  was  constructed  using  commercially 
available  radiation  software  PrismSPEC20,  calculating  the  radiation  spectrum  output  from  an  uniform  plasma  sphere. 
The  authors13  applied  the  spectra  obtained  from  plasma  sphere  with  radius  Rc  to  cylindrical  capillary  with  the  same 
radius  since  PrismSPECT  does  not  support  cylindrical  symmetry.  Although,  such  a  method  is  not  consistent 
(radiation  from  plasma  sphere  was  used  for  cylindrical  geometry),  however,  it  allowed  for  the  first  time  a  self- 
consistent  calculation  of  the  radiation  transport,  without  introducing  a  grey  factor  f  Eqs.  (1.2)  and  (1.3).  Thus,  unlike 
previous  models,  model13  did  not  use  any  “asymptotic”  radiation  models,  but  self-consistently  calculates  the 
radiation  heat  flux  at  the  transition  layer. 

The  present  paper  describes  the  high-pressure  capillary  discharge  model13  in  detail  and  extends  it  to  the  case  of 
slab  geometry  supported  by  PrismSPECT.  Analysis  of  assumptions  made  in  the  model  is  given  in  Section  II;  the 
description  of  the  model  is  presented  in  Section  III;  and  numerical  results  and  discussion  are  in  Sections  IV. 


II.  Model  Assumptions 

The  following  assumptions  are  made  in  the  model:  (1)  the  temperatures  of  electrons  and  heavy  particles  (ions 
and  neutrals)  are  equal,  i.e.,  the  LTE  condition  is  achieved;  (2)  plasma  composition  can  be  calculated  using  the  Saha 
equation,  i.e.,  local  plasma-chemistry  equilibrium  is  achieved;  (3)  heating  of  the  plasma  by  the  viscosity  drag  force 
is  small  and  can  be  omitted  from  the  model;  (4)  magnetic  pressure  is  much  smaller  than  thermal  pressure,  i.e.,  ft  < 
1;  (4)  the  the  conduction  heat  flux  is  much  smaller  than  the  radiation  flux  and  can  be  ignored  in  the  model.  The 
plasma  composition  and  geometry  assumptions  are  discussed  in  detail  in  Section  II.F. 

A.  Local  thermodynamic  equilibrium  factor  (LTE-factor) 

Since  for  typical  ablative  high-pressure  capillary  discharge,  the  mean  free  path  for  neutral-neutral  and  neutral- 
ion  collisions  are  much  smaller  than  the  capillary  radius  (capillary  slab  gap  Dc  in  slab  geometry)  the  ions  and  the 
neutrals  have  the  same  local  temperature;  they  are  in  local  thermodynamic  equilibrium. 

In  the  plasma  capillary  core  region.  Fig.  1,  where  the  plasma  is  usually  almost  fully  ionized,  or  at  least  well- 
ionized,  the  electron-ion  collisions  play  a  major  role  in  determining  local  thermodynamic  equilibrium  between 
electrons  and  heavy  particles  (ions  and  neutrals).  In  the  electron-ion  collision  process  the  colliding  particles  lose  or 
gain  some  energy  from  each  other;  the  average  exchange  energy  between  the  particles  due  to  elastic  collisions  can 
be  estimated  as 


A  ee<=>i*T 


2  •  m 


(II.  1) 
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where  T  is  the  plasma  temperature,  is  the  mass  of  an  electron,  and  M,  is  the  ion  mass.  The  average  energy  A  Wei 
that  an  electron  gains  from  the  electric  field  E  between  collisions  is 


AW.  =  ■ 


'  e,  drift 


(II.2) 


where  vei  is  the  electron-ion  collision  frequency,  and  ve, drift  is  the  electron  drift  velocity 


v  e, drift 


e-E 


m„  ■  v. 


(II- 3) 


Substituting  Eq.  (II.  3)  into  Eq.  (II.2),  we  obtain 


AW;  =  =  1.18  •  1014  •  El{V!  Wj; 1  in  [eV]  . 


m0  •  vm 


At 


74  2 

•  Z  •  n. 


(II.4) 


In  Eq.  (II.4)  we  have  used  this  expression  for  vei 


v„-  =  3.94  - 10 


TeV2[eV] 


in  [sec-1], 


(II.5) 


taken  from  the  NRL  brochure21.  Here,  Z  is  the  ion  charge,  «,  is  the  ion  density,  Acl  is  the  Coulomb  logarithm,  and  ne 
=  Z-Hi  is  the  electron  number  density  (the  plasma  is  assumed  to  be  quasi-neutral).  Now  let  us  introduce  the  LTE- 
factor,  Klte  as  the  ratio  of  Aseoi,  Eq.  (II. 1),  to  A  Wei,  Eq.  (11.4): 


1q-38  n:[m-}]-Z4-  A2 

A  Wei  '  Ai-T2[eV]-E2[V/m\ 


(II.6) 


where  we  have  substituted  T  for  Te  in  Eq.  (II. 4)  and  Mo'A-,  for  Mt  in  Eq.  (II.  1);  M0  is  the  unit  of  atomic  mass  and  H,- 
is  the  mass  of  an  ion  in  atomic  units.  When  this  factor  is  large,  the  electron-ion  collision  exchange  energy  is  larger 
than  the  average  energy  that  the  electron  gains  from  the  electric  field  between  collisions,  that  is,  local 
thermodynamic  equilibrium  is  achieved;  the  electron  temperature  is  equal  to  the  temperature  of  heavy  particles  T. 
When  Klte  <  1  the  electron-ion  collision  exchange  energy  is  smaller  than  the  average  energy  that  the  electron  gains 
from  the  electric  field  between  collisions,  that  is,  there  is  no  local  thermodynamic  equilibrium.  In  this  case  the 
electron  temperature  is  higher  than  the  temperature  of  heavy  particles  and  is  determined  by  balancing  the  energy  that 
the  electron  gains  from  the  electrical  field  with  that  transferred  to  ions.  Thus,  substituting  Te  instead  of  T  in  Eq.  (II.  1) 
and  setting  this  equation  equal  to  Eq.  (II. 4)  we  obtain  the  following  estimate  for  the  electron  temperature: 


Te[eV]  «  3  •  1CT19  • 


/?.p;r3]-Z2  ■  Aef 
A)'2  ■  E[V  !m\ 


(II. 7) 


Obviously,  the  electron  temperature  Te  obtained  this  way  forces  KLTE  to  equal  1.  In  the  case  of  multi-ion  plasma 
composition,  Eq.  (II.  6)  can  be  rewritten  as 


KLTe=  9.5-10- 


A2 


T~[eV]  -  E2[V  I m] 
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where  index  k  corresponds  to  the  ions  of  type  k.  In  this  equation  we  have  averaged  1A4,  with  the  weight  function 
equal  to  the  collision  frequency  of  electrons  with  a  given  type  of  ions. 
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and,  since  the  Coulomb  logarithms  for  all  sorts  of  ions  are  practically  the  same,  taken  Acj  out  from  the  brackets. 

Now  let  us  consider  the  LTE-factor  in  the  transition-boundary  layer,  Fig.  1,  where  the  “plasma”  temperature  is 
small,  about  0.3-1  eV,  the  gas  is  weakly  ionized  and,  therefore,  the  electron-neutral  collisions  play  the  main  role  in 
the  determination  of  the  LTE  in  this  layer.  For  the  electron-neutral  collision  process  the  average  exchange  energy 
between  neutrals  and  electrons  is  given  by  Eq.  (II.  1),  where  we  have  to  substitute  index  a  for  index  i\ 
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Here  n  and  aeak  are  the  number  density  and  the  electron-neutral  collision  cross-sections  for  type-/(  neutral  atoms. 
Substituting  into  Eqs.  (II. 2)  and  (II. 3)  the  electron-neutral  collision  frequency  veat  instead  of  electron-ion  collision 
frequency  vei  as 


v  =Yv  =VT  -Y 

ea  eak  le 
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=  6.21103  -  jTjK]  ■  Z («  [m  3 ]  •  <7  [m2]) 


where  we  have  used  an  expression  for  electron  thermal  velocity  VTe , 


VTe  = 


^8- k-T 


=  6.21  - 103  -y[rjK] 
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we  obtain  the  average  energy  AWea  that  an  electron  gains  from  the  electric  field  E  between  collisions, 
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Dividing  Aeeoa  by  A  Wea  and  introducing  the  total  pressure  in  the  transition  layer  as 

p  =  k-r-yak 

(here  we  have  neglected  the  small  electron  pressure  in  the  transition  layer),  we  obtain 


6 


(11.15) 


As  one  can  see  KLTE  in  the  transition  layers  is  independent  of  the  electron  temperature;  here  we  have  assumed  that 
the  electron-neutral  elastic  cross-sections  are  independent  of  the  electron  temperature,  a  fair  assumption  for  capillary 
discharges. 

In  our  capillary  discharge  model  we  consider  the  C4H9  polyethylene  capillary  wall  composition.  Assuming  that 
the  gas  in  the  transition  layer  is  fully  dissociated  and  using  the  electron-neutral  collision  cross-sections4 

<reC  =  2.64  •  10  19  [m2]  and  aeH  =  1.49  •  1(T19  [m2],  (11.16) 

we  obtain 

P2 

KTlte~c7h  =2.2-10  3  ■—  where  P  in  [Pa]  £in[V/m].  (11.17) 

E~ 


Thus,  if  the  LTE-factors  introduced  by  Eqs.  (II. 9)  and  (11.17)  are  much  larger  than  one,  local  thermodynamic 
equilibrium  in  the  capillary  discharge  is  achieved. 

It  should  be  stressed  that  in  the  case  of  medium  ionized  plasma  core  region,  the  electron-neutral  collisions  also 
must  be  taken  into  account  in  calculating  the  LTE-factor.  So,  Eqs.  (II.  1)  -  (II. 3)  can  be  rewritten  as 
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where  collision  frequencies  are  given  by  Eqs.  (II. 10)  and  (11.12).  Dividing  A£eoiia,  Eq.  (11.18)  by  AW,  Eq.  (11.19),  we 
obtain  an  expression  for  the  LTE-factor  for  the  case  of  moderately-ionized  plasma  core. 

In  our  model  we  calculate  the  LTE-factor,  verifying  the  applicability  of  the  LTE-approximation  used  in  the 
model.  It  should  be  stressed  that  the  LTE  approximation  was  used  in  all  previously  developed  models l  s'  1015  but, 
however,  without  clear  verification. 

II.B.  Local  plasma-chemistry  equilibrium 

Let  us  investigate  whether  local  ionization  equilibrium  in  the  plasma  core  region  of  high-pressure  capillary 
discharges  is  achieved.  Assuming  an  LTE  in  the  plasma  core  region  (the  electron  temperature  is  equal  to  the  local 
temperature  of  heavy  particles,  neutrals  and  ions,  Section  II.  A)  the  number  of  electron  impact  ionization  events  and 
the  number  of  recombination  events  in  three-body  collisions  (1  +  e  +  e  — »  n  +  e  +  h  v)  in  the  capillary  plasma  core 
region.  Fig.  1,  per  unit  of  time  can  be  estimated  as 

N7ZP  *  «.  •  Ex^~ ■  VT<e  ■  (na  ■  a,,„ )-Vc  ,  (11.21) 


Vr  recom 
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where  /  is  the  ionization  potential;  <jion  is  the  ionization  impact  cross-section  from  the  ground  electronic  state; 
(cw»a)  is  the  reverse  ionization  mean  free  path;  Vc  is  the  capillary  volume;  and  %=  8.75-10'3q-T[eV]'9  2  is  the  three- 
body  recombination  coefficient22. 

The  outflow  of  electrons  through  the  open  capillary  end  and  together  with  their  diffusion  to  the  capillary  wall 
can  be  estimated  as 


^ e. outflow  ~  ^ s 


(11.23) 


•  n  n 
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(11.24) 


where  Cs, 


C , 


r-P  _  1/ '  in,  +  na  +  ne)  ■  k  -  T 


I 0 


(n,+na)-Ma 
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is  the  sound  speed  at  the  capillary  open  end  (sonic  condition,  Fig.  1),  Ac  is  the  capillary  cross-section  area,  Rc  is  the 
capillary  radius  (in  the  case  of  slab  geometry,  one-half  times  the  slab  gap  Dc  has  to  be  used  instead  of  Rc),  AwaU  is  the 
wall  area,  Dediffus  is  the  electron  diffusion  coefficient,  Ae  is  the  electron  transport  mean  free  path,  Ma  is  the  average 
mass  of  heavy  particles,  and  T  is  the  temperature  of  plasma  core,  Fig.  1.  Assuming  the  plasma  core  to  be  fully 
ionized,  na  <  nh  Z~  \  and  substituting  y  =  5/3  we  obtain 


C 


s 


(11.26) 


When  the  rates  of  ionization  and  recombination  events,  Eq.  (23)  and  (24),  are  much  larger  than  number  of 
electrons  leaving  the  capillary  chamber  per  unit  of  time, 

N  =  N  +  N 

v  e  v  e,diffus  v  e, outflow  ’ 

we  may  conclude  that  local  ionization  equilibrium  is  achieved,  i.e.,  the  plasma  composition  is  determined  by 
balancing  the  ionization  and  recombination  processes,  resulting  in  Saha  equilibrium.  Thus,  when 


1.3  •  KT44  •  T[eVy5  ■  n]  ,  Exp\  - 


k-T 


k,  -v.o,,) 


» 


» 


rt  ni 


Avail  K 


yK  \Ma  RCJ  v 


-1.  ^  +  2.3.10'7.4M: 

’  V  M„  Rz-  A„ 


c  ee  e 


(11.27) 


the  local  ionization  equilibrium  is  achieved.  Flere  we  have  used  Eq.  (11.13)  for  VTe  and  estimated  /,L.  ~2  as 


1  «  2.3  •  10  • 
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(11.28) 


assuming  that  the  plasma  is  well  ionized,  na  ~  For  typical  capillaiy  dimensions  Rc  =  0.002  m  (Da  =  4  mm),  Lc  = 
0.1  m,  nea  =  1025  m"3,  T  <  3eV,  I/T  ~  5,  aion  ~  10"20  m2  for  carbon,  (»;/Ma)0  5  ~  0.007  (for  carbon)  and  Aee  =3we 
obtain,  that 


N“ 
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So,  for  typical  high-pressure  ablative  capillary  discharge  parameters,  the  condition  of  ionization  equilibrium  is  well 
satisfied.  Since  the  mean  free  path  for  heavy  particles  collisions  is  much  smaller  than  the  capillary  radius  for  typical 
ablative  capillary  discharge,  the  local  plasma-chemistry  equilibrium  is  also  established  when  the  local  ionization 
equilibrium  is  achieved. 

Thus,  it  has  been  shown  that  the  local  plasma-chemistry  equilibrium  approximation  is  valid  for  high-pressure 
capillary  discharges.  This  approximation  is  assumed  in  the  present  model  as  well  as  in  all  previously  developed 
models. 

It  is  worth  noting  that  in  the  transition  layer  the  photo-ionization  processes  can  be  more  important  than  the 
electron  impact  ionization  process  due  to  the  absorption  of  the  large  photon  flux  escaping  from  the  hot  plasma  core 
region.  Fig.  1,  leading  to  an  increase  in  the  plasma  temperature  in  this  layer  and  larger  plasma  density  than  that 
calculated  by  the  Saha  formula;  and  second,  the  dissociative  recombination  process  (AB ;  +  e  — >  A  +  B)  and  impact- 
three-body  recombination  process  with  participation  of  a  neutral  as  a  third  body  (i  +  e  +  a  — >  a  +  n  +  h v)  can  also 
be  more  important  than  the  i  +  e  +  e  — >  n  +  e  recombination  process  because  the  plasma  in  this  region  is  weakly 
ionized. 


II.C.  Heating  of  the  plasma  by  the  viscosity  drag  force  (VIS-factor) 

Let  us  estimate  the  plasma  heating  by  the  viscosity  drag  force.  We  can  estimate  the  thickness  of  the  viscosity¬ 
boundary  layer  using  the  momentum  equation  for  a  quasi-stationary  capillary  discharge,  Fig.  2, 
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where  p  here  is  the  mass  density  of  the  gas  in  the  viscosity-boundary  layer,  rj  is  the  viscosity  coefficient,  and  Vhoimd 
is  the  characteristic  gas  velocity  at  the  boundary  of  the  viscosity  layer,  Fig.  2. 


Viscosity-boundary  layer  V i)OUII 

Fig.  2.  Schematic  radial  distribution  of  longitudinal  plasma  velocity  in  the  capillary. 


The  heating  of  the  plasma  due  to  friction  forces  in  the  viscosity-boundary  layer  can  be  estimated  as 
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Substituting  r/~  p-  Vn'^hh  into  Eq.  (11.32),  where  VTh  and  A/,/,  are  the  averaged  theimal  velocity  and  collision  mean 
free  path  for  heavy  particles  (neutrals  and  ions)  in  the  viscosity-boundary  layer,  we  obtain 


Qvisc  ~  A wall  '  P ' 


Yjh  '  4 hh  .  y5/2 
T  "  bound 


(11.33) 


Let  us  estimate  the  rate  of  energy  leaving  the  capillary  through  the  open  end,  Fig.  1,  as 
Qoutfiow  =Ac'Cs-  Ku,  ~  A  '  Cs  ■  («,-  ■ull  +cp-k-T  ■  ( na  +  Z  • +«,.))  , 


(11.34) 


where  houi  is  enthalpy  of  the  plasma  leaving  the  capillary  per  unit  of  volume,  ud  is  the  cost  per  ionization,  cp,  is  the 
specific  plasma  heat  at  constant  pressure.  Let  us  introduce  the  viscosity  factor  (VIS-factor)  as  the  ratio  of  the 
viscosity-plasma-heating-rate,  Eq.  (11.33)  to  the  outgoing-energy-rate,  Eq.  (11.34), 
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If  KVis  factor  is  small,  then  the  heating  of  the  plasma  by  the  viscosity  drag  forces  is  small  compared  to  the  rate  of  the 
energy  leaving  the  capillary  and  can  be  neglected. 

Substituting  C,  for  Vbound,  VTi  for  VTi„  p  as  nrMa  (plasma  is  fully  ionized),  1.66  for  y,  ion-ion  mean  free  path  Au 
for  Z/,/„  and  one  for  Z,  we  obtain  an  upper  estimate  for  VIS-factor  as 
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(11.36) 


Since  the  ion-ion  mean  free  path  is  about  equal  to  the  electron-electron  mean  path,  we  can  substitute  Eq.  (11.28)  for 
An  into  Eq.  (11.36)  that  gives  for  cylindrical  capillary. 
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Substituting  typical  parameters  of  high-pressure  capillary  discharge4"7, 13:  Lc=  0.1  m,  R,  =  0.002  m,  T  ~~  3  eV,  ne  ~ 
1025  m"3,  ud  ~  15  eV,  Aee  ~  3  into  Eq.  (11.37),  we  obtain  that  KVIS  ~  0.01.  This  estimate  is  in  good  agreement  with 
numerical  results7.  In  the  case  of  slab  geometry,  we  have  to  substitute  in  Eq.  (11.37)  Dc  for  Rc. 

Thus,  it  has  been  shown  that  heating  of  the  plasma  by  viscosity  drag  forces  is  small  and  can  be  dropped  in  the 
modeling  of  high-pressure  ablative  capillary  discharges,  as  it  is  done  in  the  present  model. 


II.D.  Ratio  of  thermal  pressure  to  magnetic  pressure  (the  /^-factor) 

The  magnetic  fields  in  both  slab  and  cylindrical  geometries  can  be  estimated  as 


B\T]  = 
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and 
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where  Jsiah  is  the  total  capillary  current  per  unit  of  slab  length  along  the  y-axis,  Fig.  1,  and  /cv/  is  the  total  current  for 
cylindrical  capillary,  and  po  is  the  magnetic  free-space  permeability.  Substituting  Eqs.  (III. 3 8)  into  the  equation  for 
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we  obtain. 
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In  the  model,  we  calculate  /?  validating  a  model  assumption  of  (?«  1. 

II.E.  Ratio  of  thermal  conduction  to  radiation  heat  transfer  (COND-faetor) 

One  of  the  problems  in  zero-dimensional  modeling  of  the  capillary  discharge  is  to  accurately  estimate  the 
conduction  heat  transfer  to  the  capillary  wall  or  to  the  transition  layer  because  the  thermal  conduction  heat-transfer 
coefficient  changes  dramatically  with  the  plasma-gas  temperature,  pressure,  and  plasma  composition.  However,  in 
our  model  we  estimate  the  radial  conduction  heat  flux  at  the  outer  boundary  of  the  transition  layer.  Fig.  1,  as 


F  «  2 • 

1  T-Cond 


T-T 

Wt 

D, 


(11.40) 


where 
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are  the  electron  thermal  conductivity22  and  the  electron-electron  collision  frequency21;  T  is  the  temperature  of  the 
plasma  core,  Twau  is  the  wall  temperature,  ne  is  the  electron  density  in  the  plasma  core  region,  Aee  is  the  electron 
Coulomb  logarithm;  in  the  case  of  cylindrical  geometry  we  have  to  substitute  2-Rc  for  Dc.  Substituting  vee  into  the 
equation  for  %e  and  using  Eq.  (11.40)  we  obtain  that  total  conduction  heat  flows  at  the  transition  layer  is 
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where  the  first  equation  corresponds  to  slab  geometry  and  the  second  one  to  cylindrical  capillary.  In  Eq.  (11.42)  we 
have  dropped  Twan,  since  the  plasma  core  temperature  is  usually  much  larger  than  the  wall  temperature  or  the 
characteristic  temperature  of  the  transition-boundary  region.  Fig.  1.  It  should  be  stressed  that  Eq.  (11.42)  is  an  upper 
estimate  for  the  radial  heat  flux. 

In  the  model  we  calculate  the  ratio  of  the  radial  thermal  conduction  heat  flux,  Eq.  (11.42),  to  the  radiation  heat 
flux  at  the  transition  layer,  COND-faetor.  When  the  COND-faetor  increases  and  becomes  approximately  one,  this 
indicates  that  thermal  conduction  becomes  comparable  to  the  radiation  heat  transfer  and  has  to  be  taken  into  account 
in  modeling  capillary  discharges.  Since  in  our  model  we  assume  no  thermal  conduction,  our  model  is  valid  when 
COND-faetor  is  much  less  than  one. 

ILF.  Plasma  composition  and  geometry  assumptions 

In  our  zero-dimensional  model  of  ablative  capillary  discharge  we  make  the  following  plasma  composition  and 
geometry  assumptions: 

1.  In  the  model,  the  capillary  wall  is  polyethylene,  CrHv.  The  radiation  and  plasma  composition  database 
has  been  constructed  for  this  specific  capillary  wall  composition  and  for  4mm  slab  capillary  gap;  we 
have  used  PrismSPECT20  to  construct  our  database.  The  chosen  wall  composition  and  slab  gap  are  not 
principal;  similar  databases  can  be  constructed  using  PrismSPECT20  for  other  wall  compositions  and 
values  of  Dc. 
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2.  The  gas  in  the  plasma  core  region  is  hilly  dissociated;  this  is  valid  for  plasma  temperatures  large  than 
0.7-1  eV. 

3.  As  in  works2’8  10-15  we  neglect  the  heat  losses  into  the  capillary  wall,  assuming  that  all  radiation  heat 
reaching  the  capillary  wall  is  spent  on  its  ablation. 

4.  The  temperature  of  the  ablative  gas  entering  into  the  plasma  core  region  from  the  transition-boundary 
layer,  Fig.  1,  is  equal  to  T,  the  average  temperature  of  the  capillary  plasma.  In  other  words,  the  plasma 
radiation  ablates  the  capillary  wall  and  then  dissociates  and  ionizes  ablated  gas  and  finally  heats  it\  to 
the  plasma  temperature  in  the  thin  transition  layer,  Fig.  3.  This  assumption  has  been  made  in  all 
previous  models1"8 10-15  and,  as  it  has  been  discussed  in  Section  I,  can  be  verified  only  in  a  model  which 
would  include  radiation  heat  transfer  across  the  capillary;  unfortunately,  such  a  model  does  not  exist 
yet.  However,  experimental  observations2  4  do  not  seem  to  contradict  this  assumption. 


Fig.  3.  Schematic  of  capillary  discharge  (not  to  scale)  showing  the  transition/boundary  layer 
between  the  solid  wall  and  the  plasma.  Polyethylene  molecules  from  the  boiling 
capillary  wall  enter  the  transition  layer  (a)  and  become  dissociated  and  ionized  at  the 
boundary  between  the  transition  layer  and  the  plasma  (b).  The  temperature  of  the 
particles  leaving  the  transition  layer  for  the  plasma  core  region  is  equal  to  the 
average  temperature  of  the  plasma  core;  DBL  is  the  thickness  of  the  transition  layer. 


III.  Description  of  the  model 


III.A.  Description  of  the  database 

In  our  model  we  use  a  radiation-plasma-composition  database  that  has  been  constructed  using  PrismSPECT 
software20.  This  software  calculates  the  plasma  composition  and  the  radiative  spectral  flux  to  the  “chamber  walls” 
for  specified  total  number  density  of  atoms  (neutrals  and  ions),  the  element  composition  ratio  (in  our  case  4  to  9  for 
C4-Hg  capillary  wall  composition),  plasma  temperature,  and  geometry.  Integrating  the  spectra  over  frequency  we 
obtain  the  total  radiation  flux  to  the  walls  Frad,  which  we  use  in  our  database  as  well  as  plasma  composition 
calculated  by  the  PrismSPECT20.  Since  PrismSPECT20  offers  only  two  geometries:  spherical  and  slab  (no  cylindrical 
geometry);  we  have  selected  slab  geometry,  with  4mm  slab  gap.  PrismSPECT20  assumes  an  infinite  slab  (L,  =  oo); 
therefore,  applying  the  calculated  radiation  flux  to  the  case  of  limited  capillary  length,  Fig.  1,  is  correct,  where  Lc» 
Dc. 

Preliminary  PrismSPECT  calculations  showed  that  for  the  database  selected  temperatures  range  of  1  -  8  eV  and 
the  total  number  density  of  heavy  particles  between  1016  -  1028  m’3,  only  ///,  Hn,  Ch  Cn,  Cln,  CIV,  Cv,  Hr  (i.e.  neutral 
atoms,  ionized  hydrogen  and  singly  up  to  four  ionized  carbon  atoms)  contribute  in  plasma  composition;  the  number 
densities  of  higher  ionized  carbon  atom  are  negligibly  small  and  are  not  included  in  the  database. 
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It  is  worth  noting  that  PrismSPECT  includes  free-free,  free-bound,  and  bound-bound  radiation  and  takes  into 
account  the  pressure  and  Doppler  broadenings  of  spectral  lines  and  the  reduction  in  ionization  potential  due  to  non¬ 
ideal  plasma  effects. 

Applying  a  standard  linear-interpolation  method  to  the  database  we  calculate  the  radiation  flux  and  the  plasma 
composition  for  a  given  plasma  temperature  and  total  number  density  of  heavy  particles  (neutrals  and  ions). 

III.B.  Electrical  conductivity 

The  electrical  conductivity  plays  a  key  role  in  the  model,  by  determining  the  rate  of  Joule  heating.  Neglecting 
the  induced  magnetic  field  (/?»  1),  the  electrical  conductivity  is  given  by: 

(III.l) 

where  the  electron-neutral  momentum  transfer  collision  frequency  vea  given  by  Eqs.  (11.12)  with  the  electron 
collision  cross-sections  for  neutral  carbon  and  hydrogen  atoms  given  by  Eq.  (II.  1 6).  For  the  electron-ion  momentum 
transfer  collision  frequency  vei  we  used  the  Spitzer  equation23  modified  by  Zollweg-Liberman24  in  order  to  account 
for  non-ideal  effects. 


a  =  - 


ne 


38-Z-/7e  -e2 

rji3/ 2 

ye-me-T 


12  -n-s^-k-T 

s0-k-T 

3 

2/3\ 

/  2 

Z-e 

n  ■ e 2 

4  •  n  ■  n. 

L  J 

/ 

(III.2) 


(HE  3) 


where  So  =  8.854- 10"12  [F/m]  is  the  permittivity  of  free  space,  and  the  factor  ye( Z)  in  is  a  weak  function  of  average 
ion  charge  Z  and  can  be  approximated  as  4  ye  ~  0.58  +  0.1 -(Z-l).  It  should  be  noted  that  Eq.  (III. 2)  differs  from  Eq. 
(II.  5)  for  vei.  However,  in  the  case  where  the  non-ideal  effects  are  negligibly  small,  the  plasma  density  is  bellow  than 
1024  m'3,  Eq.  (III.2)  converges  to  Eq.  (II. 5).  Since  for  typical  high-pressure  capillary  discharges  the  no-ideal  effects 
are  rather  small  (differences  in  calculating  vei  by  Eqs.  (II.5)  and  (III. 3)  are  less  than  30%)  using  Eq.  (II.5)  to 
calculate  the  LTE-factor  is  acceptable.  However,  to  calculate  plasma  resistance  we  have  to  use  the  more  accurate 
formula  for  vei. 

III.C.  Equation  for  ablation  rate  of  the  capillary  wall,  plasma  enthalpy 

Calculating  a  required  enthalpy  hcH  to  bring  a  “ C4H9  polyethylene  molecule”  from  the  capillary  wall  to  the 
plasma  region.  Fig.  3,  is  another  key  aspect  of  the  model.  Introducing  the  ablation  rate  of  the  capillary  wall  material 
as  a  number  of  “C4H 9  molecules”  incoming  in  the  capillary  chamber  per  unit  of  time  per  unit  of  the  slab  length  along 
the  y-axis,  Fig.  1,  we  obtain 


N 


in 

CH 


2  •  F 


where  Frad  is  the  radiation  flux  at  the  boundary-transition  layer.  The  enthalpy  hcH  is 

hCH=A£„+Cp-k-T-(l  +  Z-cp)-{Ca+Ha)  , 


(HI-4) 


(HE  5) 


where  Ca  =  4,  Ha  =  9,  As#  includes  the  energies  of  vaporization,  dissociation,  ionization,  and  electronic  excitation, 
Cp  =  y!(y+  1)  is  the  specific  heat  at  constant  pressure  with  y—  5/3,  and  <p=  n/(n,-  +  na)  is  the  ionization  ratio.  The 
vaporization  and  dissociation  energies  are  considered  as  constants  in  the  model  and  are  input  parameters  in  the  code. 
It  is  worth  noting  again  that  in  the  model  we  neglect  the  radiation  heating  of  bulk  capillary  wall,  assuming  that  all 
radiation  incoming  into  the  transition  layer  is  absorbed  into  this  region  and  expended  on  the  ablation  of  wall 
material,  as  stated  in  Eq.  (III.4). 
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Calculating  electronic  excitation  energy  we  take  into  account  (1)  16,  12,  6,  4,  1,  4  and  1  pseudo-levels  (i.e. 
groups  of  elementary  levels  with  combined  statistical  degeneracy)  and  average  energy  -  for  Ch  Cu,  Cm,  Clv,  Cv,  Hh 
and  Hn  respectively;  the  contributions  of  higher  numbers  of  electron-partition  functions  into  plasma  enthalpy  is 
negligibly  small;  and  (2)  the  reduction  in  the  ionization  potentials  due  to  the  non-ideality  of  the  plasma.  In  non-ideal 
plasmas  the  ionization  potential  decreases  because  in  dense  plasmas  “free”  electrons  are  no  longer  completely  free 
(like  in  an  ideal  plasmas)  but  remain  weakly  bound  to  the  ions  (not  to  the  neutrals).  Therefore,  less  energy  is 
required  to  remove  an  electron  to  a  weakly  bond  state  (from  which  the  electron  can  then  easily  escape,  obtaining  the 
necessary  small  amount  of  energy  from  elsewhere)  than  to  a  completely  free  state.  The  reduction  in  the  ionization 
potential  of  /-ionized  atom4' 25-27  is 
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where  AD  is  the  Debye  length  and  Q  is  the  deBroglie  wavelength: 
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Here  both  the  electron  and  positive-ion  shielding  are  taken  into  account.  For  neutral  atoms,  /  is  equal  to  zero 
(they  do  not  attract  negative  particles  as  well  as  positive  ions)  so  there  is  no  reduction  in  their  ionization  potentials. 
Substituting  ne  =  6-1026  m-3,  T  =  4  eV  and  Z=  1,  into  Eqs.  (IV.  1)  -  (IV. 3),  we  obtain  AI /_>/+/  ~ /  0.3  eV.  We  would 
like  to  note  that  since  the  Coulomb  interaction  is  negative,  the  plasma  pressure  is  reduced  too  because  the  electrons 
spend  some  time  near  ions  not  bombarding  a  “chamber  wall”.  As  has  been  shown12' 27 ,  this  reduction  is  less  than 
10%  relative  to  the  total  plasma  pressure  for  our  conditions.  Therefore,  we  neglect  the  pressure  reduction  in  our 
model. 

Thus,  for  a  given  plasma  temperature  and  heavy  particle  number  density  we  find  the  plasma  composition  and 
radiation  flux  at  the  transition  layer,  Frad.  Then,  using  plasma  composition  we  calculate:  the  average  ion  charge  Z; 
the  ionization  energies  while  accounting  for  the  reductions  in  ionization  potentials  of  heavy  particles;  populations  of 
electronic  levels  of  carbon  and  hydrogen  ionized  and  neutral  atoms  by  employing  the  Saha  equation;  and  finally 
A s0  and  then  hCH.  Substituting  Frad  and  hCn  into  Eq.  (III. 3)  we  calculate  the  ablation  rate  of  wall  material. 

We  have  to  stress  that  in  Eq.  (III. 5)  we  have  neglected  the  kinetic  (not  thermal)  energy  of  the  ablated  materials 
incoming  into  the  plasma  region.  This  is  a  correct  assumption  since  the  capillary  length  L,  is  assumed  to  be  much 
larger  than  the  slab  gap  Dc,  Fig.  1. 

III.D.  Mass  and  energy  equations 

An  equation  describing  the  outflow  of  ablated  “polyethylene  molecules”  through  the  open  capillary  end,  Fig.  1, 
can  be  written  as 
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where  index  “end”  denotes  conditions  at  the  exit  plane  and  the  absence  of  the  index  “end”  denotes  conditions 
averaged  over  the  capillary  volume,  in  the  “middle”  of  capillary.  Here  iich  is  the  mole  number  density  of  ablated 
molecules  and  Mc  and  MH  are  carbon  and  hydrogen  atomic  masses  respectively. 

Combing  Eqs.  (III. 4)  and  (III. 9)  and  using  the  isentropic  flow  relations28  and  assuming  that  average  ion  charge 
and  ionization  ratio  are  preserved  through  the  capillary  up  to  the  exit  plane,  Z,,„rf  =  Z  and  (pen(i  =  (p.  we  obtain  the 
following  equation  for  the  law  of  conservation  of  mass: 
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Since  we  assume  that  there  is  no  energy  loss  in  the  capillary  discharge,  we  may  conclude  that  all  electrical 
energy  put  into  the  capillary  discharge  is  leaving  the  capillary  with  the  plasma  jet  at  the  exit  plane  or  is  expended  on 
increasing  the  inner  energy  of  the  plasma  inside  the  capillary.  This  gives  the  following  equation  for  the  law  of 
conservation  of  energy: 
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where  sCh  is  the  inner  plasma  energy, 


(III.  15) 
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Thus,  solving  Eqs.  (III. 13)  and  (III.  15)  for  a  given  Jsiab(t)  and  initial  plasma  conditions,  ncrid  =  0)  and  T(t  =  0), 
we  obtain  all  characteristics  of  capillary  discharge  vs.  time.  It  has  to  be  mentioned  that  the  authors13  used  a  similar 
system  of  equation  for  cylindrical  capillaries  investigating  the  stability  of  the  capillary  discharge  for  the  case  of 
Jsiab=cons\  and  different  initial  plasma  conditions.  They  also  calculated  the  discharge  parameters  employing  LRC 
circuit  to  self-consistently  calculate  the  current  and  compared  obtained  results  with  experiments. 


IV.  Numerical  results  and  discussions 

In  the  numerical  results  presented  below,  the  database  temperature  and  C+Hg-mole  number  density  were  1-8 
Ev,  and  7.7- 1014  -7.7T016  m-3,  the  energy  of  evaporation  and  dissociation  were  taken  as  zero  (for  polyethylene’  they 
are  much  smaller  than  ionization  potentials  of  carbon  and  hydrogen,  and  since  the  ionization  ratio  of  the  plasma  is 
larger  than  15%,  Fig.  4,  they  can  be  ignored)  The  capillary  lengths  were  chosen  as  4,  6,  and  10  cm,  and  Dc  was  4 
mm.  In  this  paper  we  investigate  the  steady-state  solutions  of  Eqs.  (III. 13)  and  (III. 15).  As  shown  in  Figs.  4  -  10,  the 
model  predicts  the  existence  of  two  steady-state  regimes  for  ablative  discharge  operation  at  a  given  plasma 
temperature.  The  first  regime  occurs  with  the  plasma  is  so  dense  that  the  radiation  mean  free  path  Ararj  is  smaller  than 
or  equal  to  about  the  capillary  gap,  the  case  of  super-high  pressure  (SHP)  capillary  discharge,  solid  curves  in  Figs.  4 
-10,  and  the  second  case  occurs  with  the  plasma  density  much  lower  such  that  Xraci  is  much  larger  than  Dc,  i.e.  the 
case  of  moderately  high  pressure  (MHP)  capillary  discharge,  dashed  curves.  As  shown  in  Figs.  4  -  10,  the  regimes 
converge  at  small  plasma  temperature,  and  there  are  no  steady-state  solutions  for  plasma  temperatures  smaller  than 
1.215,  1.43,  and  1.74  for  Lc  =  4,  6,  and  10  cm,  respectively.  In  the  SHP-  regime  the  pressure  range  is  107  —  5-  10s  Pa, 
and  in  MHP-regime  pressures  are  much  smaller  than  105  -  106  Pa,  Fig. 5,  and,  therefore  the  ionization  ratio  is  larger 
for  SHP-regime  than  for  MHP-regime  at  the  same  plasma  temperature.  Fig.  4. 

Fig.  6  shows  the  ratio  of  radiation  flux  at  the  transition  layer,  Fig.  1,  to  the  blackbody  radiation  flux  with  the 
same  plasma  temperature.  As  one  can  see  the  gray  factor f  Eq.  (1.3),  varies  from  0.1  to  almost  1  for  the  SHP-regime 
and  from  0.23  to  0.03  for  the  MHP-regime  depending  of  plasma  temperature  and  the  capillary  length.  This  indicates 
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that  the  grey  factor  may  change  noticeably  with  time  in  non-steady  operation  regime  and,  therefore,  to  assume  it  is 
constant  (as  has  been  done  in  models"'8  10  u)  can  leads  to  false  results.  The  current,  power,  and  capillary  wall 
ablation  rate  per  one  meter  of  capillary  length  in  the  y-direction  are  shown  in  Fig.  7,  8,  and  9.  As  one  can  see  the 
parameters  of  capillary  discharges  may  differ  more  than  one  order  of  magnitude  or  even  as  much  two  orders  for 
SHP-  and  MHP-regime. 


Slab  Capillary,  Dc  =  4  mm,  Lc  =  4,  6  and  10  cm 


Fig.  4.  Ionization  ratio:  —  SHP  and  —  MHP  regimes. 


Slab  Capillary,  Dc  =  4  mm  Lc  =  4,  6  and  10  cm 


Fig.  5.  Plasma  pressure:  —  SHP  and  — MHP  regimes. 


Slab  Capillary,  Dc  =  4  mm,  Lc  =  4,  6  and  10  cm 


Slab  Capillary,  Dc  =  4  mm,  Lc  =  4,  6  and  10  cm 


Temperature  [eV] 


Fig.  6.  Frad/FsBR^  —  SHP  and  — MHP  regimes. 


Fig.  7.  Discharge  current:  —  SHP  and  — MHP  regimes. 


Figs.  10  -  12  show  the  dependence  of  the  LTE-factor,  /?,  and  the  COND-factor  vs.  temperature  for  both  regimes.  As 
one  can  see  the  SHP -regime  satisfies  all  three  assumptions  of  the  model:  the  LTE-factor  and  ft  arc  much  larger  than 
one,  and  the  COND-factor  is  much  smaller  than  one  in  all  temperature  regions.  However,  the  MHP-regime  satisfies 


16 


only  two  of  the  model  assumptions:  the  LTE-factor  and  P  are  much  larger  than  one  in  all  temperature  regions,  while 
the  COND-factor  in  this  regime  increases  with  temperature  and  becomes  larger  than  1,  Fig.  12. 


Slab  Capillary,  Dc  =  4  mm,  Lc  =  4,  6  and  10  cm 


Slab  Capillary,  Dc  =  4  mm,  Lc  =  4,  6  and  10  cm 


Fig.  10.  LTE-factor:  —  SHP  and  — MHP  regimes. 


Fig.  11.  P:  —  SHP  and  — MHP  regimes. 


Slab  Capillary,  Dc  =  4  mm,  Lc  =  4,  6,  10  cm 


Thus,  the  presented  model  cannot  adequate 
describe  the  MF1P -regime  in  the  case  of  relatively  high 
temperatures,  when  the  thermal  heat  conduction 
becomes  larger  than  radiation  heat  flux  at  the  transition 
layer.  Flowever,  since  in  this  regime  the  calculated 
LTE-factor  and  P  are  both  much  larger  than  1  even 
when  COND-factor  is  approximately  1,  we  may 
conclude  even  in  the  case  when  thermal  conduction  is 
much  larger  than  the  radiation  (heat  transfer  of  the 
ablative  capillary  discharge  is  controlled  by  thermal 
conductivity  not  by  radiation)  the  capillary  discharge 
can  be  operated  at  conditions  when  the  LTE-factor  and 
(1  are  much  larger  than  1 .  This  operating  regime  may  be 
attractive  for  thruster  applications  as  well  as  SHP- 
regime  or  MTTP-regime  depending  on  specific 
applications. 

It  should  be  stressed,  that  future  investigation  of  the  stability  of  steady-state  SLIP  and  MHP  regimes  is  important, 
and  we  are  planning  to  address  this  issue  in  the  future. 


Fig.  12.  Conduction  factor:  —  SHP  and  — MHP  regimes. 
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